Microscopic model for the boson integer quantum Hall effect 
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In two dimensions strongly interacting bosons in a magnetic field can form an integer quantum 
Hall state. This state has a bulk gap, no fractional charges or topological order in the bulk but 
nevertheless has quantized Hall transport and symmetry protected edge excitations. Here we study 
a simple microscopic model for such a state in a system of two component bosons in a strong orbital 
magnetic field. We show through exact diagonalization calculations that the model supports the 
boson integer quantum Hall ground state in a range of parameters. 



Recently there has been considerable theoretical ac- 
tivity on the possibility of "Symmetry Protected Topo- 
logical (SPT)" insulating phases of interacting bosons in 
diverse dimensions. These are generalizations of the cel- 
ebrated free fermion topological insulatorspQ to systems 
of bosons. As in the non-interacting limit a boson system 
is trivial, strong interactions are necessary to reach any 
insulating phase, including the topological ones. Thus 
studies of bosonic insulators force us to think of topo- 
logical insulation without the crutch of free fermions and 
band topology appropriate for non-interaction fermion 
systems. It is expected that the insights from these stud- 
ies will be useful in thinking about the interplay of strong 
correlations with topological insulation in electronic sys- 
tems. 

In the interacting context we seek phases of bosonic 
matter that have a bulk gap, have no 'intrinsic' topo- 
logical order or fractionalized excitations, but neverthe- 
less have protected boundary excitations. In space di- 
mension d = 1 the familiar Haldane spin chain pro- 
vides a well established example of such a state. Re- 
cently progress was made in classifying all such gapped 
lei systems based on the concept of group cohomology[5]- 
5 . In d > 1 a cohomology classification has also been 
proposed[5] of all gapped bosonic systems with no 'in- 
trinsic' topological order which predicts the existence of 
many SPT states. Physical properties of these phases 
are not readily accessed through this framework and have 
been addressed through other methods in both twofTHS] 
and three cfimensions 10 KT|. States beyond the coho- 
mology classification are also known [TU1 IT21 1X1] to exist 
(ind = 3). 

One of the most important questions that remains 
open in this field is the identification of experimentally 
relevant microscopic models for these SPT phases in 
d > 1. Here we address this question for a prototyp- 
ical SPT phase in d = 2, namely an integer quantum 
Hall state of bosons. We show that a simple and realis- 
tic Hamiltonian for a system of two-component bosons in 
a strong magnetic field provides a concrete realization of 
this state. The physical properties of this state have been 
studied in Refs.mUnH andQH The electrical Hall con- 



ductivity a xy was shownjS] to be quantized to be an even 
integer (see Ref. |9] for a simple argument) . At the edge 
to the vacuum there is a charged chiral edge mode and 
a counter propagating neutral mode. Despite being non- 
chiral this edge structure is protected so long as charge 
conservation symmetry is preserved. This boson inte- 
ger quantum Hall state provides a prototypical example 
of a Symmetry Protected Topological insulator in boson 
systems in two space dimensions. Ref. [§] proposed that 
two component bosons in strong orbital magnetic fields 
at filling factor v = 1 for each component provide a suit- 
able platform to realize the boson integer quantum Hall 
state (bIQH). A number of potential competing states 
were also identified. 

In this paper we study the simplest Hamiltonian for 
such a system through exact diagonalization studies. We 
find strong evidence for the stability of the boson integer 
quantum Hall state. 

We consider a bilayer system of N interacting bosons 
and N$ flux quanta. All our calculations are performed 
in the lowest Landau level and the filling factor is defined 
as v t = N/N$. The bilayer index is equivalent to a spin 
1/2. The number of particles per layer is respectively iV^ 
and N^. We define the spin projection as S z = (N' — 
N^)/2. We denote by zj,zf the complex coordinates of 
the two boson species. The Hamiltonian projected in the 
lowest Landau level reads: 

h = 9s j:s^ (4-4)+^ (2) (4-4) 
+<? d £* (2) (4-4) a) 

When g s — g c i the Hamiltonian has an extra pseudospin 
SU(2) symmetry which rotates the two species of bosons 
into one another. The ground state of this model was 
previously studied using exact diagonalization in Refs. fTTl 
and [TH] with an emphasis on total filling factor vt = |. 
Here we focus instead on v t = 2 which is a suitable plat- 
form to realize the bIQH state. Ref. |9] showed that pro- 
totypical wave functions for the bIQH state is a pseudo 
spin singlet. Thus we may expect the bIQH state to 
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be stabilized, if at all, near the g s — gd point. A com- 
peting pseudo spin singlet state is a member of a fam- 
ily of Non-Abelian Spin Singlet (NASS) statespS] pro- 
posed for two-component bosons at various total fillings 



2k 



v t 

the bIQH state wins over the NASS state at total filling 
factor Vt = 1 + 1 = 2. When the ratio 2^ = 0, the two 
components are decoupled and at v = 1 each. With delta 



function repulsion there is strong evidence [20 22 that 
each component is in the non-Abelian Pfaffian (Moore- 
Read) state (Pf)E3]. This decoupled Pf x Pf state of the 
2-component system will be stable in some range of small 
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The corresponding wave function reads 
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oo we expect phase separa- 



In the other limit when iLsL 
tion into puddles with each puddle consisting of bosons 
living in one of the two layers with twice the average 
density (hence at filling factor v — 2). Thus we expect 
the region around — = 1 to be most favorable for the 
bIQH state and to see phase transitions to Pf x Pf or 
the phase separated state as this ratio is varied. Our 
calculations below confirm this expectation. We will see 
that there is a sizable gap at the SU(2) symmetric point 
which persists in a wide range of the ratio — . Thus the 
bIQH state seems quite robust in this model The phase 
diagram obtained from our work is summarized in Fig. 
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FIG. 1. Phase diagram of the model Hamiltonian. The bIQH 
state is stable in a range of — estimated to be between ~ 0.8 
and ~ 1.3 from exact diagonalization calculations with 12 and 
16 particles. PS refers to the phase separated state which 
appears for large — . For small — the decoupled Pfaffian 
state Pf x Pf is stable. 



A wave function for the bIQH state was proposed in 
Ref.l 



^ flux — P, 
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where Plll denotes the projection onto the lowest Lan- 
dau level. As there is no intrinsic topological order there 



is a unique ground state on the torus. 

The competing NASS state can be built up starting 
with the (to, to, n)-Halperin wave function given by: 



(k integer). Our calculations below show that ^ 
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The filling factor is v t = 2/ (to + n). Note that the 
(1, 1, 1)-Halperin is the densest zero energy state for our 
model Hamiltonian in Eq. [TJwhen g s ,gd > 0. From the 
Halperin wave function, one can build a series of NASS 
states at filling v t — 2 m-i ■ This is done by dividing the 
particles into k groups, writing a Halperin (to; to, to — 1) 
state for each group and then symmetrizing over the dif- 
ferent groups. This procedure leads to: 

fc-i 

*m = <S[II +1> • • ' > Z (i+1) g ' 

i=0 



(5) 



where S is the symmetrization operator. For m = 2, n — - 
1, these are the NASS states introduced by Ref. [19]. 
The NASS state is ( fc + 2 H fc + 1 ) degenerate on the torus 
geometry. 

For our numerical calculations, we consider the torus 
geometry with a square aspect ratio. The choice of the 
torus instead of the sphere allows us to avoid the bias of 
the so-called shift: In many cases, such a bias prevents 
the direct comparison of possible candidates in finite size 
calculations. Moreover, each model state has a well de- 
fined non-trivial degeneracy on the torus, providing a 
simple signature that we can extract form the energy 
spectrum. We performed calculations using the transla- 
tion symmetry of the torus along the y direction. Thus 
our states are only labeled by the k y momentum which 
can be expressed as an integer between and iV$ — 1. We 
first consider the SU(2) symmetric point (i.e. g s = gd). 
We have looked at several system sizes up to 18 parti- 
cles. A typical energy spectrum is shown in Fig. [2]a, for 
N = 16 at filling factor v = 2. We clearly observe a 
unique ground state at k y — and with total spin S = 0. 
Such a signature is the one expected for an IQH state at 
this filling factor. There is no trace left of the Pf x Pf 
that should occurs in the limit gd = 0. Its 9-fold degen- 
erate manifold (5 states in momentum sector k y = and 
4 states in momentum sector k y = 4) is completely lifted. 
All the energy spectra at v = 2 that we have computed 
share exactly the same features as this N = 16 example. 
In Fig. [3j we give the behavior of the neutral gap A (i.e. 
without changing the number of particles or the number 
of flux quanta) as a function of the system size. While 
a thermodynamical extrapolation is not straightforward 
due to the size effects, it looks convincing that A will 
stay finite. Note that contrary to most of the other FQH 
phase, the gap seems to slightly increase with increasing 
N. 
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The hallmark of the FQHE is the unique nature of 
its excitations which are realized by inserting or remov- 
ing flux quanta or particles. In the case of the integer 
quantum Hall at v = 1, removing on electron leads to a 
number of degenerate ground states that exactly matches 
the number of orbitals. For many of the model states 
of the FQHE, allows to predict how many states should 
appear in the low energy spectrum using the Haldane's 
exclusion principle [21] • Each counting is a signature of 
the considered phase. Fig. [2]d shows the low energy part 
of the spectrum when we remove one particle from the 
case described in Fig. [2^,. There is a set of low energy 
states with spin S = ¥, one per momentum sector. In 
this example, we thus get 16 states. This is exactly what 
would be obtained for the bIQH state when removing 1 
particle from a finite system 16 fcrmions. Similar results 
are obtained when adding one particle. This is a clear 
indication that the low energy physics of our model is 
equivalent to the bIQH state. Other ways to nucleate 
excitations such as adding or removing flux quanta, gen- 
erate much more states. Due to the energy splitting in 
finite size calculations, the signatures in the energy spec- 
trum are not as clear as in the case of adding removing 
one particle. 

We have tried to look at the particle entanglement 
spectrum[551 [55] of the ground state. Unfortunately, 
we were not able to extract any clear signature. Note 
that the bIQH phase is a composite fermion-like state 
with reverse flux attachment and such a behavior is ex- 
pected: Such states (like the fermionic spin-polarized 2 /3 
state) have a particle entanglement spectrum that looks 
a generic wave function for the same system size (i.e. the 
corresponding reduced density matrix does have neither 
an exponentially large number of zero eigenvalues nor an 
entanglement gap). We have also looked at the sphere 
geometry at the shift given by the model wave function 
of Eq. [3] For every system size, we have observed a spin- 
single state with a total angular momentum equal to zero 
and a large gap. 

We now turn to the phase diagram when tuning the 
ratio gd/gs- Figs.|4^, and|4]D for two system sizes, N = 12 
and N — 16. Note that the Pf x Pf phase occurs when 
N is a multiple of 4, and the NASS state occurs when 
N is a multiple of 6. So N = 12 is the only case where 
we can have the three phases. In the gd/g s — limit, we 
observe a 9-fold low energy manifold that corresponds to 
Pf x Pf phase. However, one cannot distinguish it for 
N = 16. There is finite splitting of the 3-fold degenerate 
MR ground state when this phase is realized through the 
two-body delta interaction. Thus when having two copies 
of this system in finite size, the splitting of the total 9- 
fold low energy manifold can be larger than the gap of a 
single copy. 

The NASS state with k = 3 and m = 2 is a potential 
candidate at v t = 2. A first signature of this state should 
be a 10-fold low energy manifold on the torus (the degen- 
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FIG. 2. Left panel: Low energy spectrum for N — 16 bosons 
and N<s> — 8 flux quanta on the torus geometry. The ground 
state is a spin singlet (S = 0) clearly separated from the 
higher energy excitations. Right panel: Low energy spectrum 
for N = 15 bosons and N& = 8 flux quanta on the torus 
geometry. This situation corresponds to one less particle that 
the left panel. Since N is odd, the spin is half integer. The 
exact degeneracy is a consequence of the magnetic translation 
symmetry when N and N<s> are relative prime numbers. 
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FIG. 3. Neutral gap A as a function of 1/N at filling factor 
v = 2. While a thermodynamical extrapolation is quite diffi- 
cult, the trend strongly suggests that the gap will stay finite 
for large N. 



eracy being only approximate since we are not using the 
model interaction that produces this state). For N = 12, 
we have never found such a signature for the full range 
of gd/g s values. This should be contrasted with v t = | 
where Refs. [T7I and IT81 found old evidence for the k = 2 
NASS state. 

In the phase diagram, the most prominent phase is 
related to the one discussed for the SU(2) symmetric 
point gd/g s — 1- For both system sizes, the gap is sta- 
ble over a large region of gd/g s that include this spe- 
cific point. To get a more quantitative estimate, we 



can compute the overlap 
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I 2 between 



the ground state any value of gd/gs with the ground 
state at gd/gs = 1- For TV = 12, this overlap is higher 
than 0.9 for 0.80 < gd/g s < 1-25 (i.e. centered around 
the peak in Fig. |4^l). For N — 16, this range becomes 
0-75 < gd/gs < 1-40 (for this system size, the Hilbert 
space dimension is 5.2 x 10 6 ). At the peak value of gap 
around gd/g s — 1-7, the overlap is still 0.78. 

In summary we have studied a simple microscopic bo- 
son model that shows strong evidence for an integer quan- 
tum Hall ground state with no intrinsic topological order. 
This state is a prototypical example of a Symmetry Pro- 
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FIG. 4. Phase diagram for TV = 12 (upper panel) and N = 16 
(lower panel) on the torus geometry at filling v — 2 when 
tuning the interaction parameter gd/g s - The red line indi- 
cates the gap AbiQH between the lowest energy state and the 
first excited state. Each kink in AbiQH corresponds to a level 
crossing above the ground state. The blue line stands for 
the gap Apfxpf between the low energy manifold of 9 states 
associated to the Pf x Pf phase and the first excited state. 
In both cases, the gap is set to zero if low energy manifold 
that corresponds to each of these states, does not have the 
proper quantum numbers. The shaded area indicates the re- 
gion where the lowest energy state is fully polarized (i.e. a 
phase separated state). 



tected Topological insulator of bosons in two space di- 
mensions, and is an interacting boson analog of the free 
fermion topological insulators. The existence of this in- 
teracting dominated topological insulator phase in such 
a simple model provides hope that other such bosons in- 
sulators in both two (6] [8] and three [10] dimensions may 
also be realized in fairly simple models. 

Clearly a natural experimental context to seek a re- 
alization of the boson integer quantum Hall state is in 
ultra-cold atoms. The delta function repulsion is realis- 
tic and natural in cold atom systems. Realizing an or- 
bital magnetic field big enough to be in the quantum Hall 
regime is a long standing challenge in the field. We hope 
that the results of this paper provide further impetus to 
meet this challenge. 

Note: During the preparation of this letter, we became 
aware of two related papers[27l [28]. In particular, we 
recover the same scaling of the gap on the torus geometry 
than Ref. |27]. 
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